Entanglement of low-energy excitations in Conformal Field Theory 
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In a quantum critical chain, the scaling regime of the energy and momentum of the ground state 
and low lying excitations are described by conformal field theory (CFT). The same holds true for 
the von Neumann and Renyi entropies of the ground state, which display a universal logarithmic 
behaviour depending on the central charge. In this letter we generalize this result to those excited 
states of the chain that correspond to primary fields in CFT. It is shown that the n-th Renyi entropy 
is related to a 2n-point correlator of primary fields. We verify this statement for the critical XX 
and XXZ chains. This result uncovers a new link between quantum information theory and CFT. 
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Entanglement is one of the central concepts in quan- 
tum physics since Schroedinger used the term in an an- 
swer to the Einstein-Podolsky-Rosen article in 1935. A 
particularly active line of research is concerned with the 
role played by entanglement in the physics of many-body 
systems [ij. One is typically interested in the amount of 
entanglement between two spatial partitions, say A and 
B, of a many-body system in its ground state. For a 
pure ground state the amount of entanglement is usu- 
ally quantified with the entanglement entropy, or the 
von Neumann entropy of the reduced density matrix pA- 
Sa = — trA PA In PA- Alternatively, the Renyi entropies 
Sn are also used: {Sn)A = j^lntr^ Pa", the entan- 
glement entropy being lim„^i Sn- One of the most im- 
ortant results in this topic is the celebrated area law 0^ 
which, roughly speaking, states that ground states of 
gapped many-body systems with short-range interactions 
have an entanglement entropy proportional to the area 
of the hypersurface separating both partitions. The area 
law restricts the fraction of the Hilbert space accessible to 
ground states of local Hamiltonians in an essential way, 
allowing for their efficient numerical simulation 

Violations of the area law occur in gapless (critical) 
systems. In one dimension most of critical systems, as 
well as being gapless, are also conformal invariant. The 
attention to the entanglement properties on these sys- 
tems came after the seminal result of Holzhey, Larsen and 
Wilczek , who showed that the leading behavior of the 
ground state entropies S^^ is proportional to the central 
charge of the underlying conformal field theory (CFT) 
governing the long-distance physics of the discrete quan- 
tum chain. If £ and N are the lengths of the partition A 
and of the total system, both measured in lattice spacing 
units, then the Renyi entropy of _the_ ground state, with 
periodic boundary conditions, is 
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where c is the central charge of the CFT and 7„ is a non- 
universal constant. 

In a critical model, the finite-size scaling of the en- 



ergy of excitations is given by the scaling dimension of 
the corresponding conformal operators Q . This fact sug- 
gests that also the entanglement entropy could be related 
to properties of these operators. Entanglement of excited 
states has been considered previously. In [§] it was shown 
that the negativity of the excited states in the XXZ criti- 
cal model shows a universal scaling. In (lo| it was shown 
that a violation of area law should be expected for the 
low lying excited states of critical quantum chains, and 
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, it was considered the entanglement of very large 
energy excitations in XY and XXZ spin chains. 

In this letter we show that the entropy 5'°'^'^ of excited 
states associated to primary fields exhibits a universal 
behaviour that generalizes (dJ . The energy of these low- 
lying states degenerate as 1/iV in the bulk limit TV — ?> oo. 
We prove that the excess of entanglement, S^'^ — S^^, is 
a finite-size scaling function related to the 2n-point cor- 
relator of the primary field. These results are verified in 
two models: the XX and XXZ spin chains. 

Entanglement of generic primary states. Let us 
consider a system S of length N with periodic bound- 
ary conditions. To describe it, we introduce the complex 
variable C = cr + it, where < o" < iV is the spatial 
coordinate and t is the time coordinate. S is split into 
two subsystems S ~ AU B, with A = {e,£ — e) and 
B = {£ + e,N - e), and where e « £ < N is a. short- 
distance cutoff 0]. The world sheet of the past (t < 0), 
is a cylinder with two semidisks of radius e cut out (de- 
noted C and D in figure [T]) . The boundary of the world 
sheet of figure [T] is given by the union AU C U B U D. 
After the conformal transformations C — >• w — >■ z: 
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the C cylinder gets mapped into a strip of height tt and 
width d = 21og [;^sin(^)]; being A, B, C, D the 
boundaries of the strip in z space (see figure [T]). More- 
over, the point at the infinite past Coo — —ioo gets 
mapped into Coo — !• Zoo = in {1 — £/N). We shall con- 
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FIG. 1. Riemann surfaces describing the past events in and 
z. The distinguished point in z is the infinite past (^oo = — ioo. 



sidcr the simplest excited states in a CFT, namely the 
primary states, which are those generated acting on the 
vacuum |0) with a primary field T(C, C), with conformal 
weights {h, h), 



|T)= _lim T(C,C)|0). 



(3) 



The wave function of this state |T) is given by the path 
integral 



*xy(T) (X V(j)T 



(4) 



where (p denotes the local field whose Euclidean action is 
S {(/)). The field T is a functional of </>, that is evaluated at 
the infinite past Zoo in equation (|4]) (recall equation (O). 
X and Y denote the values of the field </> in the subsys- 
tems A and B respectively. Periodic boundary conditions 
are imposed on the C and D edges ■ If T were not pri- 
mary, then equation ^ would include additional terms 
generated by the conformal transformations ([2]). 

The density matrix p = pA of subsystem A is obtained 
by tracing over the variables in B: 



Pxx'(T)cx JvY *xy(T)*^x,(T). 
Plugging (j4]) into (O one finds 



(5) 
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where z'^ ~ i7r(l + £/N) represents the point at the 
infinite future. The functional integral is over a strip 
of height 2tt and width d, with boundary conditions 
(j) = X on the lower edge and (j) ~ X' on the upper edge. 
The normalization factor is determined by the condition 
trp = 1, which implies that Z(l) is the functional inte- 
gral with no operator insertion and the top and bottom 
edges of the strip being identified (i.e. a torus partition 
function), and (TT^) is the two point correlator on this 



torus. To compute the entanglement entropy one first 
)mputes the trace of p^, which is given by 
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where Z[n) denotes the partition function on a torus of 
lengths 27rn and d, so that the moduli parameter is given 
by T„ = 27rm/c?, and where (. . .)t-„ denotes the expec- 
tation value in the r„-torus. Notice that the 2n-point 
correlator of fields T, depends on the ratio t/N and 
on the moduli parameter. 

To further proceed one uses the expression of the par- 
tition function Z{n) of a general CFT with central 
charge c for chiral and antiehiral sectors of the theory, 
Z{n) = Z(T,r) ~ tr (7^f'~24 q^'>~2i^ with the nome 
q = q = exp (27rzT). In the limit d >> 1, it is convenient 
to perform the modular transformation t — > —^/t. The 
partition function is modular invariant and can be easily 
evaluated in terms of the nome q = e^'^'^^l'^ = "^Z". In 
particular, for the ground state (Tq ~ 1) one gets (up to 
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as anticipated in ([1]). In the general case, equation ([7]) 
depends on a 2n-point correlator of the fields T and T^^ 
on a cylinder of length 2?™ along the time direction. It is 
now convenient to rescale this length to 27r. Afterwards, 
we shift the coordinates Zj — ^ Zj — i7r(l — x)/n where 
X = £/N. Finally, we exchange a and t coordinates in 
such a way that zj ~ li^j jn for T and Zj ~ 2TT{j + x)/n 
for T^. The ratio between the excited and the ground 
state traces, F^\x) = tr p^/tr py^, becomes, from ([7]): 
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where (. . .)cy denotes the expectation value in a cylinder 
of length 2tx. Note that F^^ ^ exp [(1 - ?7,)(S'J - S^^)] . 
The dependence of the entropies of the excited states 
on the n-point correlation functions was also observed in 
the ground state entropies of two disjoint segments of the 
quantum critical chains [l2j . The entanglement entropy 
for the excited state |T) can then be computed using the 
replica trick: 
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In the hmit x << 1. the terms TT''' appearing in ^ 
can be approximated by the operator product expansion 
(OPE) TxTt = l + 1' + ... , finding: 
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{x)^l + ^ (--n) {nxf + 0{x'^- ) (11) 



where is the operator with the smahest scahng di- 
mension, Aij,. The term of order a;^^* depends on 
the OPE constants C*^! and on the expectation val- 
ues («'(0)*(^))cyi. If A* = 1, this term is 0{x^) as 
the first one in equation pip , and eventually they may 
cancel one another, as we shall see in an example below. 
If A,!, ^ 1 one could use pT|) to infer the quantities h + h, 
A^jj and C^^-^ from the numerical computation of the en- 
tanglement. 

Using equation (|10l) one finds, for the low- a; behaviour 
of the entanglement entropy {£/N << 1): 



2 f t ^ 2A4, 



(12) 



Equations (|WT^ are the main results of this letter. 
They relate the von Neumann and n- Renyi entropy of the 
excitation represented by the primary operator T to the 
2ri,-point correlators of T and in the cylinder. Notice 

(n) 

that the ratio F-^ does not depend on the non-universal 
constant 7„, which is therefore common to and S'^'^ . 

As an example of the laws (|9ll0p we shall consider a 
c = 1 CFT given by a massless boson compactified on a 
circle. The primary fields are given by the vertex oper- 
ators Ti[n,m] = e*("+''^+"-'^) (being 0,0 chiral and an- 
tichiral boson fields) where a± = n/2^/K±my/K, k is the 
compactification ratio, and n, to G Z. The scaling dimen- 
sions of these operators are (a^ -|-a?_)/2 = n'^ /AK + rn^K. 
Using the chiral correlator of vertex operators on the 
cylinder, ([{, e^''^^(^-^^)cy = n,>fc P sin(z,fc/2)]-"^"^ it 
turns out that 



Ti[i,fe] 



(x) = 1, yx,j,k. 



(13) 



Hence, all the excitations represented by vertex operators 
have the same entropy as the ground state. This result 
is not in contradiction with (|lip because, in this case, 
A^ = Agcj, = 1 and both 0{x^) terms in ([TT|) cancel out 
due to the properties of the OPE constants. In fact, the 
cancellation happens in all order of x. 

Next, we study the operator T2 = idcj). Us- 
ing its correlator on the cylinder (90(zi)90(z2))cy = 
— [2 sin(zi2/2)]~^ and the Wick theorem, we get (in terms 
of s{x) = sin(7ra;/2)): 



and a more lengthy expression for i^^"^^' • ^^e \aw-£/N 
limit, one finds that Fy'\x) '--^ 1 -I- (ttx)^ (l/n — n) /3, 
which leads to an excess of entanglement entropy given 
by dni) with {h,h) = (1,0). 

Realizations of both types of excitations in particular 
models will be now shown, and their amount of entan- 
glement compared with the CFT predictions (|13ll4p . 
Excitations in the XX and XXZ models. The 
Hamiltonian of the spin-1/2 XXZ model is given by 



1 ^ 



), (15) 



where N is even and periodic boundary conditions are 
assumed (for A = we get the XX model). This model 
is integrable [l3| and gapless for — 1 < A < 1. The cor- 
responding CFT is given by the aforementioned bosonic 
CFT with K = f [7r-cos-i(-A)]"\ The XX model 
in the sector with magnetization M ~ ^ erf can be 
mapped, through a Jordan- Wigner transformation, into 
a system with np — M + N/2 free fermions in a lattice 
of N sites. We computed the entanglement and Renyi 
entropies of several types of excitations in these models. 
This task was achieved using the methods of references 
0, in the free fermion problem and through numerical 
exact diagonalization in the XXZ case. 

Let us first consider the vertex operator Ti[0,to]. 
In the free fermion model, the result (|13p is ex- 
act and can be proved analytically. Indeed, 
|Ti[0, to]) corresponds to the umklapp excitation 

njli4f+(2j-i)VJv'^-'=^-+(2j-i)Vw|0)' ^here kp = 
•nnp/N is the Fermi momentum, and where |0) is the 
Fermi state and d\ the fermionic creation operator with 
momentum k. This state can be obtained from the Fermi 
state shifting all the momenta as fc — ?> A: -I- 27x11: /N . Such 
a shift produces a global phase factor in the wavefunc- 
tion in real space and, consequently, the entropy remains 
unchanged. In the XXZ model, the state |Ti[0, 1]) cor- 
responds to the ground state in the sector with np spins 
up and total momentum P — 2'Knp/N . We observe that 
the prediction [x]^ = 1 holds, up to the oscillations 



expected for n > 2 [IJ] , which in this case are of the or- 
der of 10"'^ for systems with = 30 spins. 
We will now consider the excitation Ti [2, 1]. In a system 
of free fermions the resulting state corresponds to the ad- 
dition of two fermions at the right of the Fermi point, i.e., 

|0). Figure [2] shows that 



to the state dt ,_/^r^i^ 



Fi^l\x) = 1 - 2s^{x) + Ss-^ix) - 2s^{x) + s^{x) (14) 



kF+Tr/N^kF+3TT/Ni 

ground and excited states entropies 5*2,3 coincide, up to 
oscillations. In the XXZ model, |Ti[2,l]) is the low- 
est eigenstate with total momentum P ~ 2Ti{np -f 2)/N . 
Again in this case, oscillations of F~^/^r^ around one are 
observed (see figure [2]). 

Finally, figure [3] shows some numerical results for the 
entanglement of the excitation T2 = idcf). In the free 
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X = i/N 

FIG. 2. An illustration of the law F^^^2 ~ 1]^^^ ~ 

for the XX model with N = 100, nj? =^50, and for the XXZ with 
A = -1/2 and AT = 30, np = 14 (16 for the excited state). The 
entropy of ground and excited states coincide and follows the law 
|[T]| (continuous lines), up to oscillations [TS ). 
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FIG. 3. The quantity for the XX model at different filling 

fractions and for the XXZ (A = -1/2, N = 30, np = 14) model, 
vs. the CFT prediction l|14|l . For np = 250 the oscillations around 
I I14II arc so small that both curves are indistinguishable. In the 
inset we show S'2 for the ground and excited states {np = 250, 
N = 500) . The upper inset is a zoom of the region selected by the 
small rectangle over the curve in the main figure. 



fermion problem, IT2) corresponds to a particle- hole ex- 
citation: dJ^^_^^^j^dkp^TT/N\0) , while in the XXZ model 
it corresponds to the lowest eigenstate with P = 2tt/N. 
We observe an excellent agreement with the theoretical 
prediction ([M)) for n — 2. Similar results hold for n — 3. 
Moreover, we have checked, for n up to 6, that the low- 
£/N formula ([T^ is very well satisfied for fermions. 

In summary, we have obtained an expression for the 
Renyi entropies of excitations associated to any primary 
field. We verified the results with finite-size realizations 
of the XX and XXZ models up to 30 sites in the latter 
case, finding very good agreement with the theory. 



As explained earlier, equation ^ can be used as a nu- 
merical method to extract information about correlators, 
conformal dimensions and OPE coefiicients of primary 
fields. An interesting problem is to generalize these re- 
sults to the descendent states in CFT. We expect that the 
Renyi entropies, at a given level of a conformal tower will 
depend on the particular state targeted. This can pro- 
vide a method to establish the correspondence between 
degenerated excited states of a critical lattice model, and 
the descendent fields in the underlying CFT. 

Equation ^ further suggests a generalization of the 
Renyi entropies in terms of traces of different density 
matrices tr [pxi/OXa • ■ •]• This object would be related to 
the correlator: (T1TIT2T2 . . .) in the very same fashion 
as in (jni). The numerical computation of the associated 
generalized entropies would then provide information on 
more general correlators in CFT, and vice-versa. Ap- 
plications of the present work to other models and to 
non-primary fields are in progress. 

This work represents a further step along the direction 
of deriving CFT data using quantum information meth- 
ods. 
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